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Abstract 

It is proposed that the prompt emission observed in bursts that 
exhibit a thermal component originates from relativistic radiation me- 
diated shocks that form below the photosphere of the GRB outflow. 
It is argue that such shocks are expected to form in luminous bursts 
via collisions of shells that propagate with moderate Lorentz factors 
r < 500. Faster shells will collide above the photosphere to form col- 
lisionless shocks. We demonstrate that in events like GRB 090902B 
a substantial fraction of the explosion energy is dissipated below the 
photosphere, in a region of moderate optical depth r < 300, whereas 
in GRB 080916C the major fraction of the energy dissipates above 
the photosphere. We show that under conditions anticipated in many 
GRBs, such relativistic radiation mediated shocks convect enough ra- 
diation upstream to render photon production in the shock transition 
negligible, unlike the case of shock breakout in supernovae. The re- 
sulting spectrum, as measured in the shock frame, has a relatively low 
thermal peak, followed by a broad, nonthermal component extending 
up to the KN limit. 



1 Introduction 



According to the standard view, the prompt GRB emission is produced be- 
hind intern al shocks that form in the coasting region of a baryon loaded 



fireball (e.g. iLevinson &: Eichlerl . Il993l ; iRees fc Meszarosl . 11994 ; ISari fc Piran 



19971 ). The emission mechanism commonly invoked is synchrotron cooling 



of relativistic electrons that are accelerated at the shock front, perhaps ac- 
companied by synchrotron-self Compton emission at highest energies. In 
order to be able to accelerate particles the shocks must be collisionless. As 
argued below (see also Levinson & Bromberg 2008, hereafter LB08), under 
the conditions anticipated such shocks can only form in regions where the 
optical depth for Thomson scattering is smaller than unity, as otherwise the 
radiation produced in the shock transition and in the immediate downstream 
mediates the shock via Compton scattering. Thus, an implicit assumption 
of the standard fireball model is that a significant fraction of the bulk energy 
dissipates above the photosphere. 

The fireball model outlined above was originally motivated by the detec- 
tion of nonthermal spectra in many GRBs. Despite some difficulties in asso- 
ciating the low energy part of the spectrum with optically thin synchrotron 
emission (e.g. Crider et al. 1997; Preece et al. 1998), so me believed that 



19931 ) inferred 



this model can account for the Band spectrum (IBand et al 
in many sources. 

The difficulties with the synchrotron model, as well as other considera- 
tions (e.g. Eichler & Levinson 2000), have led to the suggestion that photo- 
spheric emission may contribute to the prompt GRB spectrum, in addition to 
the observed non-thermal component (e.g. Meszaros & Rees 2000, Meszaros 
et al. 2002). At the same time, it has been argued that a combination of a 
thermal and a nonthermal components can fit better the data in quite a num- 
ber of bursts (e.g. Ryde 2004, 2005, 2006, Pe'er et al. 2007, Pe'er 2008). The 
recent detections of a prominent thermal peaks in severa l LAT bursts, no- 
tably GRB 090902B flAbdo et all l2009bl : lRvde et allboioh . provide a strong 
support to this interpretation. 

The conditions at the shock formation region depend on the power and 
Lorentz factor of the outflow and the duty cycle of the central engine. For 
shells to collide above the photosphere, the Lorentz factor must be large 
enough. As shown below, in some bursts, e.g., GRB 080916C, internal shocks 
are expected to form above the photospher, whereas in others, e.g., GRB 
090902B, a major fraction of the energy is likely to dissipate below the pho- 
tosphere, albeit in a region of moderate optical depth. I ndeed, the spectrum 



of GRB 080916C is well fitted by a broken power law flAbdo et all 12009a) 



while that of GR B 090902B exhibits a thermal peak flAbdo et all . [2009b 



Rvde et all 120101 ) . 



The above considerations motivate detailed investigation of the properties 
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of shocks that form below the photosphere, in regions where the Thomson 
optical depth exceeds unity. Preliminary studies indicate that such shocks 
are mediated by Compton scattering of radiation prod uced in the shock 



transition layer or c onvected by the upstream flow (LB08; iKatz et all 12010 



Budnik et all l2010l ). The typical width of the shock transition, a few Thom- 



son mean free paths, is much larger than any kinetic scale involved, rendering 
particle acceleration in such shocks highly unlikely. Thus, the photon spec- 
trum produced in relativistic radiation mediated shocks (RRMS) may be 
considerably different than that anticipated in collisionless shocks. More- 
over, the radiation produced downstream is trapped for times much longer 
than the shock crossing time, and diffuses out after the shock breaks out of 
the photosphere and becomes collisionless. Consequently, any rapid fluctua- 
tions of the flow parameters will be smeared out, and may not be imprinted 
in the observed gamma-ray emission. 

In this paper we explore the role of RRMS in the prompt phase of GRBs. 
A qualitative discussion of RMS is given in §2. In §3 we analyze regimes 
in the parameter space in which overtaking collisions of shells occur below 
and above the photosphere. In §4 we investigate the properties of RRMS 
that form during the prompt GRB phase. We show that under conditions 
anticipated in many GRBs, such relativistic radiation mediated shocks con- 
vect enough radiation from the upstream to render photon production in the 
shock transition negligible, unlike the case of shock breakout in supernovae 
and hypernovae (e.g., Weaver 1977; Katz et al. 2010). Bulk Comptonization 
then produces a relatively low thermal peak, followed by a broad, nonthermal 
component in the immediate downstream. To illustrate some properties of 
the spectrum we present, in §5, test particle Monte Carlo simulations of bulk 
Comptonization. The shock structure computed in LB08 is used as input to 
these simulations. In the regime where photon convection is important the 
analysis of LB08 is justified, and the shock solutions obtained there present 
reasonable approximations. We find that a hard tail extending up to the KN 
limit (roughly ~ Tm e c 2 in the observer frame) is a generic feature of RRMS 
in GRBs. We conclude in 56. 
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RMS - a qualitative discussion and review 
of previous work 



The st ructure o f non-relativistic RMS w as first studied by lPail (119661 ) . IZerdovich &: Raizer 
(119671 ) and by IColgate Chen! (119721 ). These authors assumed that the 
plasma and the radiation field maintain local thermodynamic equilibrium 
(LTE) at the shock transition, giving rise to a black body spectrum. Under 
this assumption the tempera t ure in creases monotonically from the upstream 
to the downstream. IWeaverl (119761 ) extended this analysis to include condi- 
tions anticipated in supernovae shock breakout, where the upstream is cold 
and photon poor. He showed that when the shock velocity exceeds ~ 0.03c 
and the upstream temperature KT U < 1 n\^ b eV, where n^m is the baryons 
mass in units of 10 15 gr, the photon production rate is insufficient to sustain 
a black body distribution within the shock transition. As a consequence, the 
photons can only maintain Wien equilibrium with the plasma, defined as a 
Planck spectrum with a non-vanishing chemical potential, and a tempera- 
ture equal to that of the electrons. Since fewer photons supply the shock 
pressure, the temperature rises above its value at the far downstream, where 
thermodynamic equilibrium is re-established. At (3 U ~ 0.2 the peak tempera- 
ture reaches values of KT ~ 10 KeV, approximately two orders of magnitude 
above the black-body temperature further downstream. At such velocities 
the photon production rate in the transition region is negligible, and the 
majority of photons which support the shock are creat ed downstream o f 
the deceleration zone, and diffuse back to the upstream ( jKatz et all 120101 ) . 
Higher upstream velocities lead to increased pair production- and annihila- 
tion rates up to a point where the amount of MeV photons begins to affect 
the flow. In this case, the full KN cross section m ust be used to corre ctly 
solve the shock. RMS in this limit are analyzed by ( IBudnik et all 120101 ) . 

When the photon-to-electron ratio in the upstream region is high enough, 
the total production rate of photons in the transition region (due to local 
emission processes and diffusion from the downstream) is negligible; it may 
be assumed that the photon number is conserved within the shock. Blandford 
& Payne (1981a,b; hereafter BP81a,b) investigated such shocks in the limit 
of non relativistic upstreams and vanishing thermal effects, and obtained a 
self-consistent solution for the structure of the shock and for the transmitted 
radiation spectrum. They showed that photons which are advected with 
the flow across the shock, can diffuse back to the upstream and traverse 
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the shock once more. In each crossing the photons gain a small amount of 
energy as a result of the flow compression. Some of these photons are caught 
in the shock, and by experiencing this process multiple times, undergo a 
substantial energy buildup, resulting in a power law tail which extends to 
high energie s, similar to first order Fer mi acceleration of cosmic rays (for 



a review see iBlandford &: Eichlerl . 119871 ) . The solution provided by BP81b 
for the transmitted spectrum is valid as long as the average energy gain per 
scattering due to compressional heating (bulk motion), ~ (3 2 , sufficiently 
exceeds the corresponding energy gain due to thermal Comptonization ~ 
-^j. This condition is equivalent to the requirement that the Compton 
parameter satisfies yg < 1, and holds for fast, c old flow s . BP 81b s olution 
was g eneralized by iLyubarsky fc Sunyaevl (119821 ) , iRiffertl (119881 ) and iBecker 
(119881 ) to include thermal processes as well. 

In the non-relativist ic regime, where the upstream velocity /3 U <C 1 and 
the temperature across the shock T <C m e c 2 /K, several approximations can 
be made which greatly reduce the complexity of the solution. Firstly, the 
Thomson cross section ot can be used in the radiation transfer equation, 
instead of the full KN formula. Secondly, as the radiation field anisotropy 
in the fluid rest frame is always small, the motion of the photons across the 
shock can be regarded as a diffusive process. In such process the number 
of scatterings that a photon undergoes as it passes through the shock decel- 
eration region is N = (L s /Xt) 2 , where L s is the width of the deceleration 
region and At = (n e or) -1 is the photon mean free path. If the diffusion time 
roughly equals to the shock crossing time then \ T N ~ L s j f3 u . Consequently 
iV ~ f3~ 2 and L s ~ (n e o"r/3 u ) _1 . Finally, to order (3 2 , the radiation field 
satisfies the equation of state P r = U r /3, which provides a closure condition 
for the set of hydrody namic equations governing the shock structure (e.g. 



Hsieh & Spiegell . Il976h . 



When (3 U approaches unity the photon distribution function becomes 
highly anisotropic and the diffusion approximation breaks down. A com- 
pletely different method of solution is then required. In addition, pair cre- 
ation inside the shock transition may be important and needs a proper treat- 
ment. LB08 computed the structure of the shock transition for various up- 
stream condition, using multiple moments expansion of the transfer equation. 
They assumed that the photon number is roughly conserved across the shock 
transition, which is justified when photon advection by the upstream fluid 
dominates over photon production, and neglected pair creation by nonther- 
mal photons. They found, as expected, that the thickness of the shock is 
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of the order of a few Thomson lengths. Recently, iBudnik et al.l (120101 ) cal- 
culated numerically the structure and spectrum of RRMS for sufficiently 
relativistic, cold upstream. They have shown that the spectrum consists of 
a thermal peak at m e c 2 roughly, followed by a hard tail that extends up 
to r^m e c 2 , as measured in the shock frame. The high temperature in the 
immediate downstream is a consequence of the paucity of photons advected 
from the upstream, resulting in intensive pair cascade at the immediate down- 
stream which keeps the temperature at this level. A much lower temperature 
is anticipated in the RRMS which are expected to form in the photons rich 
jets of GRBs, as discussed below. 



3 Formation of RRMS in GRB outflows 

To investigate the conditions under which shells may collide below the pho- 
tosphere we consider a conical fireball having an isotropic equivalent lumi- 
nosity Li SO , We assume that the fireball is ejected with an initial Lorentz 
factor To ~ 1 from a compact central engine of radius Ro, and that it carries 
baryons with an isotropic mass loss rate M&. The properties of the fireball, 
and in particular the location of the photosphere depend of the dimensionless 
parameter rj = jf^- When r] < rj c , where 

1/4 

\4ir R m b c 3 

the fireball is sufficiently opaque, such that the radiation is trapped during the 
entire acceleration phase. The major fraction of the explosion energy is then 
converted into bulk kinetic energy of the baryons, and the fireball reaches 
a terminal Lorentz factor ~ 77 at some radius r* ~ r]R /T , beyond 
which it continues to coast. The photosphere is located somewhere in the 
coasting region, at r p h > r*. On the other hand, when r\ > r\ c the fireball will 
become transparent already during the acceleration phase, before reaching 
the coasting radius r = 7] c Rq/Tq. The Loerentz factor in that case is limited 
by i] c 0. 



^ = ( ) = 1-8 x 10%{^X f \ < 1 > 



1 If 77 > rj c the acceleration continues beyond the photosphere, since the photon flux is 
high enough to sustain efficient acceleration even though only a fr action of the photons 
undergo scatterings (see Meszaros fc Reest bOOOt iNakar et al~ . 2005 ). 
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The optical depth at a radius r, defined as r(r 
expressed in terms of the fireball parameters as 



f r °° <7 T n[T 1 dr, can be 



r(r) 



vtRo f [3 



2(vJv)T 1 (vR /r r) 2 Vc<V, 

1 V < Ve- 



il) 



The photospheric radius r p h can be found from the condition T(r p h) = 1, and 
it is readily seen that for r\ = r\ c one has r ph = r] c R /T = r*. 

Suppose now that intermittencies of the central engine lead to ejections 
of shells that collide at some radius r^. Let Ti denote the Lorentz factor of 
a shell ejected at time to and = bTi the Lorentz factor of a second shell 
ejected at t + St. If b > 1 the shells collide at r d ~ 2T\c-J^-r5t. Now, if 



rj < r] c then the shells collide in the coasting region and we have Ti 
optical depth at the radius of collision is obtained from Eq. ([2]): 



7]. 



T[Td) 



2T oV 5 



c5t 

Rq 



^ fl 



The 



(3) 



and it is seen that collision will occur below the photosphere, viz., r(rd) > 1 
if 



Rq 



(4) 



As can be seen from Eq. (j4j) it is easier for internal shocks to form below 
the photosphere in bursts with high luminosity. Such high luminous bursts 
with L iso ~ a few 10 53 ergs/s are detected now by the Fermi satellite. These 
bursts are characterized by an initial prompt phase with a maximal energy 
< 10 MeV and a peak energy at about 1 MeV, followed by a longer phase 
characterized by a harder spectrum extending up to 10 GeV and a larger 
peak energy. The Lorentz factors, inferred from the requirement that the 
shell should be optically thin for the GeV phot ons, are usually high with 



T > 500, but these are model dependent (see e.g. IZou et all 1201 ll . for lower 
values). Table 1 shows some relevant values for 3 such bursts, where T min i is 
the minimal estimated Lorentz factor as suming two sepa rate emission zones 
for the MeV and for the GeV photons (IZou et all l201ll ). and T min ^ is the 



minimal Lorentz factor assuming a single emission zone. 



2 For simplicity we assumed an infinite medium. If the shoc k forms inside a finite 
shell, the upper limit of the integration s hould change to 2r (e.g. IShemi fc Piranl [199(1 



Meszaros et all Il993t iNakar et all 120051 ). but that will have only minor effect on our 
results. 
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In fig (pQ) we plot r(r<i) as a function of the shell's Lorentz factor, for the 
3 bursts shown in table 1. The bold lines depict constant 5t, and span a 
range between Rq/c and a maximal time interval, taken to be the observed 
variability time by the Fermi BGO detector. The colored areas mark the 
range of T and St relevant for each burst. As can be seen shells with T > 700 
collide above the photosphere, resulting in shocks which are most likely col- 
lisionless. On the other hand, slower shells ejected over time scales which 
are closer to the dynamical time scale of the system, collide below the pho- 
tosphere leading to the formation of RMS. The parameter space explored in 
fig (0Q) indicates that in GRB 080916C, the majority of the GRB energy dis- 
sipates above the photosphere, therefore a thermal component may be very 
week or absent. In GRB 090902B, however a major fraction of the energy 
is likely to dissipate below the photosphere, in a region of moderate optical 
depth (r < 300). Indeed a prominent thermal component in the prompt soft 



phase accompanied by a hard tail was reported in this burst (lAbdo et al. 



2009bh . as expected in a case of RRMS (see In the short GRB 090510, 



a moderate fraction of the jet energy may dissipate below the photosphere. 
The interpretation of the soft initial phase in this burst however is not clear 
and may fit a Band spectrum such as expected from Collisionless shocks, as 
well as a photospheric emission with a (quasi) thermal component (see e.g. 
Ackermann et al. . 2010 ; Pelassa fc Ohno . l2010l ). 



Note that much longer variability time scale is anticipated in a case where 
the shock forms below the photosphere, and the downstream radiation is 
trapped for times much longer than the shock crossing time, and diffuses 
out after the shock breaks out of the photosphere. Therefore the observed 
variability timescales do not necessarily reflect the actual time separation 
between the shells, which may be much shorter. At what depth downstream 
of the shock equilibrium is established is yet an open issue. Since the enthalpy 
downstream is dominated by radiation, a full equilibrium is not expected for 
the moderate optical depths found above (see further discussion in §5]). The 
photons produced by RMS can reach energy of > Tm e c 2 = 0.25 (5^5) GeV 
in the observer frame, which is more then enough to account for the energy 
range observed during the soft phase. 
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Table 1 Parameters of Fermi LAT bursts 



GRB 


^90 (s) 




>rg) 


Liso 


(erg/s) 


Kar W) 


L min,l 


L min,2 


080916C 1 


66 


8.8 • 


10 54 


3.3- 


1Q 53 


512 


488 2 


890 1 


090510 3 


0.5 


1.1 ■ 


1Q 53 


2.7- 


1Q 53 


14 


324 2 


720 3 


0909025 4 


30 


3.63 


• 10 54 


2.4- 


1Q 53 


53 


253 2 


550 5 



(1) lAbdo et all (l2009ah: (2) IZou et all dpOllh : (3) lAckermann et alJ (|2010h : (4) 
Abdo et alJ (|2009bl l: f5 Ulvde et all <|20ld > 

* The variability is measured with the BGO detector (150 Kev - 40 MeV). 



4 Properties of RMS with photon rich up- 
streams 

Consider an infinite plane parallel RMS moving in the x direction. The shock 
transition is defined within a region bounded by x s — L s and x s , where = (3 U 
if x < x s — L s and /3 = (3d if % > x s . Far from the transition region the pho- 
tons maintain thermodynamic equilibrium with the plasma, and have a black 
body distribution with an average energy per photon (hv') u ^ = 3KT^ d y 
Hereafter primed quantities are measured in the fluid rest frame and unless 
stated otherwise, unprimed quantities refer to the shock frame. We define 
the shock upstream to be photon-rich, if the "effective" increase in the pho- 
ton number when passing through the shock transition is small. The term 
" effective" is used to indicate that only those photons which signi ficantly con- 



tribu te to the pressure are accounted for (see similar analysis in iKatz et al. 
201(1 but in a different regime of parameters spance). Taking T' s = T'(x s ) 
as the typical temperature at the shock transition^!, photons which are pro- 
duced with energies <C KT' S can also contribute to the radiation pressure, if 
they can be IC upscattered to energies of ~ KT' S during the passage through 
the shock. In non-relativistic RMS, where thermal Comptonization is the 

2 

dominant heating mechanism, a photon can double its energy after ~ ™kt' 
scatterings. Therefore, the lowest frequency above which the photon energy 



3 In shocks where bulk acceleration is significant, the spectrum is non-thermal; defined 
as such, this temperature represents the ratio between the radiation pressure and the 
photon number density. Note that the peak temperature in the shock may be higher than 
T' s , yet we find that in the parameter regime relevant here, the difference never exceeds an 
order of a few. 
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2 

can significantly increase is set by the c ondition «££Arj|§g < 1, where a// 



is the free- free absorption coefficient (e.g. IRybicki &: Lightmanl . Il979l ). giving 



hu' min ~ 50 < 5 



KT' a 
lOOeV 



-5/4 



Such a photon can reach the peak energy 



within the shock crossing time if yc > ln[ , , 8 ) (IWeaveij, Il976l ; iKatz et al. 

V min / 

2010). If, however, bulk acceleration is the dominant mechanism, yc — 1, 
and only photons with energies above hv' min ~ KT' s /2 can contribute to 
the pressure. When 7 U > 1 each photon undergoes only a few scatterings, 
experiencing an energy increase < 7 2 . Therefore /M^ n ~ KT' s j^\ at this 
limit. 

In addition to the Comptonization processes, emission can also be an im- 
portant source of ph otons if the phot ons production time is shorter than the 
shock crossing time. ISvenssonl (119841 ) demonstrated that the two most domi- 
nant emission processes in a thermal plasma, when pair density is negligible, 
are e-p Bremsstrahlung and double Compton (DC). The effective photon 
generation rate by Maxwellian e-p Bremsstrahlung emission is: 

-1/2 

(njer T c)7i£aA//, (5) 



■ / 

n 



ff 



32vr / V 



m P cr 



where A 



-Ei 



f) ) • 9ff, Ei(x) is the exponential integral of x, and 



gff is the Gaunt factor (e.g. IRybicki fc Lightmanl . 119791 ) . For the above es- 
timations of hv' min it is safe to assume that A// < 10 for most upstream 
conditions. Here and for the rest of the paper, electrons and protons 
sumed to have a Maxwellian distribution. This assumption relies on the fact 
that the relaxation time of an electron in restoring a Maxwellian distribu- 
tion with temperature T e is of the orde r of t e ~ (tijGt c lnA) _1 (kT e /m e c 2 ) 3 / 2 , 
where A is the Coulomb logarithm (see ISpitzerl . Il978l . BP81a). For electrons 
with subrelativistic temperatures this time is much shorter tha n the aver- 



age time between scatterings, t c ~ (n^a-rc) -1 . ISvenssonl ( 11984J ) calculated 
the production rate of photons by DC emission in the limit of soft photons 
[hv 1 min (m e c 2 , KT')), having a Wien distribution and maintaining ther- 
mal equilibrium with electrons. He showed this rate to be 



h' 



DC 



3tt 



(/iV 2 > 



m P cr 



n'A 



DC, 



hi 



where n' r is the photon number density, A DC - ^> . , ... 
is a numerical factor which replaces the gaunt factor. We stress that the 



KT' 



(6) 

g DC , and g DC 
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photons within th e transition region may not be in local eq uilibrium with 
the electrons (e.g. Blandford fc Paynd . Il981bl ; iRiffertl . Il988h : nevertheless, 



Chluba et al.l ( 120071 ) showed that if the average photon energy exceeds that 
of the electrons, the emission efficiency is reduced. Consequently, eq. ([6]) can 
be used as an upper limit. 

Last, diffusion from the downstream provides an additional source of 
photons to the shock transition region. Downstream of x s the photon number 
density continues to rise due to emission processes, until the photons regain 
thermodynamic equilibrium with the plasma. This leads to propagation of 
photons backward toward the upstream in a manner which can be treated as 
a diffusive process as long as h' r /n' r « (cft^ d (3d,c). The distance, L D , which 
photons can cover by diffusion against the flow before being swept back, can 
be estimated by maximizing the expression Ld = Id — Pdd, where Id = 



2yDt is the diffusion length and D = 3n ? aT is the diffusion coefficient (e.g. 



Weaver! . 1 19761 . BP81a). This gives Ld = 3n , \ T p d i implying that photons can 



diffuse back as long as (3d < 1/3, which always holds in a shock downstream. 
The total number density of photons at point thus be estimated by 

summing the total number of photons swept from the upstream with those 
which diffuse back from the downstream. 



„, <^ L + Ct° s !^ (1 + Q , (7 ) 

Ud \ K,u U u J U d 

where h' r = hjf + n' DC . The dimensionless parameter ( quantifies the rela- 
tive contribution of photon production inside the shock and at the immediate 
downstream. Specifically, £ -C 1 implies that the net production rate is much 
smaller than the incoming flux of photons, so that the radiation field in the 
shock transition is dominated by photons advected from the upstream. The 
upstream in this case is defined to be photon-rich. Relativistic temperatures 
or the presence of sufficient amount of pairs introduce additional emission 
processes such as pair creation and annihilation and e + e_ Bremsstrahlung. 
However, as shown below, our analysis is restricted to situations where pairs 
are unimportant and temperatures are well below m e c 2 , therefore this possi- 
bility is ignored. 

An upper limit on £ can be derived analytically using the following as- 
sumptions: 

1. Within the shock transition the photon flux is conserved to a good 
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approximation, i.e. n' r u u u = n' rd Ud and 

KT ' S = Pr,d u d/K jU u u . (8) 

2. The photon production rate at the shock transition and the immediate 
downstream is constant and equals to the peak production rate of the 
shock at x s . 

3. The photons at x s have a Wien distribution with a temperature that 
is equal to that of the electrons. In this limit (/i 2 i/ 2 ) = 12 (KT' S ) 2 . 

The upper limit of ( resulting from these assumptions is: 

C c = (o.o2A^' M er 3 er°- 5 A // + 7&^A DC ) , (9) 

where O' = -^-C. As it is shown below each term in the brackets becomes 
important at different regimes of upstream conditions. The first term is 
generally more important in the non relativistic limit, while the second term 
becomes important when the upstream has relativistic velocities. 

In the limit where emission processes are negligible, the photon spectrum 
in the shock transition and in the immediate downstream is determined by 
a combination of thermal and bulk Comptonization. When thermal Comp- 
tonization is the dominant process, the photons establish a kinetic equilib- 
rium with the electrons, resulting in a Wien spectrum. A hard non-thermal 
tail may evolve only when photons are upscattered by the bulk motion of 
the fluid at the upstream, and cross the shock before being thermalized. 
For a population of electrons with a non-relativistic temperature, the pho- 

ton thermalization time is tx ~ ( n b arc J j . The shock crossing time 
is estimated differently for non-relativistic and relativistic upstreams, where 
in the non-relativistic case t s ~ (3n' 6 o"rc/ /3 2 ) ■ If however the upstream is 
relativistic t s ~ (^aTcr/7 2 ) , where r is the optical depth of the shock 
transition for a photon coming from the downstream. LB08 calculated r to 
be of the order of a few (r ~ 3 — 5), implying that a mutual parameter 
for both cases can be defined, T = ~ 12 — 1^5. If T < 1, a photon 
that is upscattered at the upstream can cross the shock without loosing its 
energy to thermal scattering, and a high energy tail may evolve. Numerical 
simulations we preformed show that in the limit of T — > 0, this tail is hard 
and reaches energies up to 7 u m e c 2 in the shock frame (see $5]). 
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4.1 Limits on the upstream conditions 

In order to analyze the conditions under which the shock upstream can be 
considered photon-rich, we express the downstream quantities in eq.flH]) in 
terms of the upstream values. Those are derived from the fluid equations, 
assuming conservation of total mass, energy and momentum across the shock. 
When /3 U <C 1 the shock jump conditions can be written as: 

n'b,uPu = n' btd p d , (10) 

(n' b , u m b c 2 ^ + \p' r ^ /3 U = {n' b4 m p c^ + ^ fa, (11) 

K,um p c 2 (3 2 u + P r ; u = n' b4 m p c 2 (3 2 d + P^ d , (12) 

Here we assume that at the upstream and downstream the contribution of 
baryon pressure is negligible. Such assumption is easily justified in photon- 
rich upstreams, since the number density of photons far exceeds that of the 
electrons and protons. In terms of the upstream parameters we therefore get: 

^H i+8 wk)' (i3) 

PU-\<,urn^-\^^- (14) 

Note that the radiation pressure in the upstream is limited by P' ru < |n' b u m p c 2 f3^. 

Above this limit (3 U is below the soundspeed, and no shock can form. We 

pi 

now substitute P' r u and n' b u with the dimensionless parameters P 



r,u 



b,U r 



,,m v c 



and N = -FS which relate to the original upstream parameters through: 



b,u 



KT' U ~ 90 eV P- 3 N-\ (15) 
' b u ~ 2.4 x 10 15 cm- 3 Pl 3 N^ 4 . (16) 



The resulting shock characteristic temperature in the non-relativistic limit 
is: 

KT>, NR = 0.2 KeV jfr^ l + igg^ (17) 
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and consequently, the derived upper limit for the relative number of newly 
created photons within the shock is: 



(_ 2 x -2.5 
1 + 0.8P_3/Ti j +12x 10 -4^- 2/3 2 iAdci _ 

(18) 

Since both parameters are defined in terms of the upstream values we dropped 
the subscript u for convenience. 

When 7 U /3 U > 1, T' s and ( c can be estimated by solving the relativistic 
generalization of eqn. fllOH12j) and taking {3d ~ 1/3. This gives: 

KT' s R ~ 23 KeV A^i (} + 0.04P_ 2 ) , (19) 

and 

C R ~ 0.1iV- 2 7l 3 A DCil 1 + 0.04P_ 2 1 + 1 + 0.04P_ 5 

v / 7l" A DC,1 v 

(20) 

The numerically calculated values of T s ' and £ c as a function of 7 U /3 U , 
A" and P are plotted in Fig (j2J). The figure is divided into 6 panels which 
represent different cuts in the N — P plane. The left hand panels show 
plots made with P = 10~ 3 , 10 _1 , 10. In the right hand panels we fixed TV 
to be iV = 10 3 , 10 4 , 10 5 from top to bottom respectively. The orange lines 

/ jap/ \ 

are contours of equal T s in units of log 10 ( -^f- J , whereas the light blue lines 

represent contours of equal ( c , with values ( c = 0.1, 1, 10. Below the black 
line the upstream Mach number < 1, and no shock occurs. The region 
where photon advection from the upstream prevails in the shock transition 
is highlighted in gray, and is bounded by ( c = 1 and KT' S < 0.1m e c 2 . At this 
temperature the energetic photons in the tail of a Wien distributed radiation 
constitute < 1% of the total photons number, and their effect may be ignored. 
The real energy distribution may also include a non thermal energetic tail, 
where is such KT' S represents the average photon energy rather than 

the Wien peak. This may affect the result, especially in the high 7„/3 u regime 
(see $5]). The green line describes constant values of /ogio(T), where values 
< imply that the photons don't have time to reach thermal equilibrium with 
the electrons in the shock transition, and a high energy tail evolves. This 
allows for the generation of a high energy tail in the immediate downstream. 



It can be seen that there is a considerable regime of parameters that result in 
RMS with photon-rich upstreams which are expected to develop prominent 
non-thermal tails. 



4.2 Applications to GRB jets 



Shocks may form in GRB jets due to various reasons, such as collision between 
converging flows (e.g. iLevinson fc Eichlerl . l2000t iBromberg fc Levinsonl . 120071 ) ; 
interaction of t he jet head with a stellar envelope of a collapsar, prior to its 
breakout (e.g. iMatznerl . l2003t IBromberg et all 1201 lh : or collisions of con- 
secutive relativistic shells. But no matter the source of the dissipation, if it 
occurs below the photosphere it will form a RMS. Here we use the results 
obtained above to show that these shocks will most likely have photon-rich 
upstreams, and generate photon spectra with relatively low thermal peaks. 
Specifically in the case of shell collisions, the upstream velocity is mildly 
relativistic and the spectrum is expected to have a hard non-thermal tail. 

Suppose that a shock is formed at some optical depth, r, below the pho- 
tosphere of a fireball, such as the one discussed in §3. The shock upstream 
condition is determined by the properties of the fireball at the correspond- 



ing dissipation radius. We can therefore use eq. 
parameters, N and P at that radius giving: 



i.25 x lo^^-r- 1/4 , 



Vc 



to calculate the flow 



(21) 



Pr 



T 



1/3 



4/3 



3-2 2£ 



1/3 



r l/3 , ^ 



4/3 



V < Vc- 



(22) 



By placing these expressions in eqn. (jTT]) - (|20p we can extract limits on the 
values of r\ which render the shock upstream photon-rich at that optical 
depth. In case where the shock has relativistic upstream, we get: 



660L^- 1/2 7«,ir; /2 



<„<i2 X i() 4 L 1/A n l/2 'y- 3 r 1/2 T- 1 

V ~ L - Z X iU -^52 n 6 !uA L T ■ 



(23) 



Below and above this regime the temperature at the immediate downstream 
exceeds ~ 50 KeV which we set as the limit above which and pairs may affect 
the flow. In the non relativistic limit shocks can form only if P < 0.75/3^, 
corresponding to an upper limit of 



V < 480L^*< / V 



-1/4 



(24) 
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where the lower limit should be derived numerically. Fig. ([3]) depicts the 
numerically calculated region for ( c < 1 in the ^ U (3 U vs. t]/t] c plane. Here we 
put r = 1 and Tq — 1 for illustration. The color codes are similar to fig.(T5]). 

In the standard internal shocks model where shocks from from shell colli- 
sions, t and 1] are related via Eq. ([3]). The requirement that the shocks form 
below the photosphere restricts the jet's baryon load to rj < r) p h, defined in 
eq. (j4|). This sets limits on N and P, where taking 10 2 7] 2 < r] < rjp^ gives: 

e.g-io^VVo 172 < N grb < 2.5.10 4 f^^V /5 ^ 2 1/20 ^ /10 r 1/2 , 



(25) 



b 2 c5t 
and 

31 ' 10 - 3 (^^) 8/15 ^ 2/15 ^ /I5 < p ~ < ^ (^f> 

(26) 

where b = T2/T1 is the ratio between the Lorenz factors of the colliding 
shells, and Ti = 77. If we further assume that the shells have roughly the 
same mass, we can parameterize the upstream Lorenz factor at the shock 
frame. Equal mass shells have a center of mass which moves with a Lorenz 
factor Tcm — ^ivb. Consequently the Lorenz factor of each shell relative 
to the central mass is Ti :C m = ^i^cm (1 — PiPcm) = ^t|> which for b > 1 
gives an upstream Lorentz factor at the shock frame 

In < V2T 1<CM = i±^. (27) 
v2o 



(2? 



Placing these constraints in eq. (120]) gives 

2/5 



6 • R- 3/5 T b^ (A) < Cc < OMLl^VTob^ 



which imply that the shock transition is dominated by the upstream photons, 
and its temperature, derived from eq.( !T9|) is: 

6.5KeV Lli 20 R^ w VT-b ( S^) ^ < KT' S < 23KeV L^i^V 



\ Ro ) 

(29) 

At such a temperature T < 1, and therefore a hard non-thermal tail may 
grow beyond the thermal peak. Such radiation-mediated internal shocks are 
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therefore expected to exhibit an observed thermal peak at energies of the 
order of a few MeV, together with a prominent hard non-thermal tail which 
may extend to a few hundreds MeV, depending on the optical depth at the 
shock location. A similar spectrum was observed in the initial phase of GRB 
090902B. 

5 Bulk Comptonization: Test particle MC 
simulations 

To illustrate some basic properties of the spectral energy distribution of 
RRMSs with photon rich upstreams, we have preformed a 2D test particle 
Monte Carlo simulations, using the shock profiles computed self-consistently 
by LB08. For convenience, an analytic fit to these profiles was adopted, 

U(0 = 0.5 [tanh«) - 1] [U d - U u \ + U d . (30) 

where U u is the far upstream 4-velocity, Ud is the far downstream 4-velocity 
and a is a free parameter that determines the width of the transition layer. 
The dimensionless variable £ is related to the angle averaged optical depth 
across the shock. In terms of the coordinate x it is expressed as 

£ = n u a T x, (31) 

where n u is proper density far upstream and cr T is the Thomson cross sec- 
tion. As stated above, these global shock solutions ignores KN effects and 
pair productions, and are suitable in the regime where convection of seed 
photons by the upstream flow is sufficiently large. As will be shown below, 
the neglect of pair production may not be justified at all for large enough 
upstream Lorentz factors. We have kept the full KN cross section in our test 
particle MC simulations to examine its effect on the spectrum, but ignored 
the thermal spread of the electrons in the shock and downstream. Thus, the 
results are relevant for photons above the thermal peak. 

The injected photons are drawn from a monochromatic and isotropic (in 
the fluid rest frame) source of photons, located at a fixed position within 
the simulation box, as viewed in the shock frame. We have run simula- 
tions for different source locations, including injection from the immediate 
upstream and immediate downstream, and found little differences in the re- 
sultant spectrum as long as the optical depth from the injection point to the 
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shock transition is not too large. The boundaries of the simulation box are 
defined by two values of the variable £, £, min and £, m ax, where £ min is located 
upstream and £ max downstream. The boundaries £, m i n and £ max are chosen 
to optimize the run time without affecting the resultant spectrum. That is, 
it has to be large enough to allow injected photons to undergo the maximum 
number of scattering across the shock, but minimum number of scattering in 
the far downstream, where the probability to recross the shock becomes very 
small. The location of £ m j„ is chosen such that less then 2% of the injected 
photons escape the simulation box though the upstream boundary. 

For each of the results presented in Figs. (TjJ and (jSJ) below 10 7 photons 
were injected in the immediate downstream. The trajectory of each photon 
was followed until it "escaped" the simulation box from either side. Fig. @ 
delineates the effect of the upstream Lorentz factor on the spectrum at the 
immediate downstream. The width of the shocks in both cases is r = 3, 
which roughly gives the profile computed in LB08, where r is the Thomson 
depth for a photon moving from the downstream, at U(£d) = Ud(l + 10~ 4 ), 
to the point U(£ u ) = 0.9U U upstream, viz., 



with 7(f) = [f/ 2 (£) + 1] 1/2 and f3 = U/^y for the profile given in Eq. EDJ 
As seen, the spectrum hardens with increasing Lorentz factor, but even for 
modest ones we find that a considerable fraction of the shock energy is trans- 
ferred to photons near the KN limit. The effect of the shock width on the 
photon spectrum is explored in Fig. |5j We have also performed runs for 
different injection points, including photon injection upstream, and found 
only small differences in the resulted spectra. We conclude that the location 
where photons are injected does not affect the spectrum considerably. 

In all cases studied a hard spectrum was obtained, extending over several 
decades in energy. For sufficiently high injection energies the statistics al- 
lowed us to follow the spectrum up to the KN limit, •y u m e c 2 , as measured in 
the shock frame. There is no evidence for a cutoff or steepening at energies 
well below the KN limit, in any of the experiments performed. As naively 
expected, the spectrum is somewhat harder for larger upstream Lorentz fac- 
tors. There is also a weak dependence on the shock width. The profiles 
computed in LB08 are relatively narrow, in the sense that the angle aver- 
aged optical depth of the order of a few (r ~ 3 — 5). As demonstrated by 




(32) 
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Budnik et al.l ( 120101 ) . pair production upstream can lead to a photon bread- 



ing cycle that can boost the upper cutoff well above the KN limit, up to an 
energy of 7^m e c 2 in case of sufficiently relativistic shocks. The temperature 
in the immediate downstream in the case of photon r ich upstream is muc h 
smaller than in the case of a cold upstream studied by iBudnik et al.l (120 10[ l. 
Nevertheless our test particle simulations illustrate that pair production by 
the accelerated photons is expected to be important, especially in the case 
of moderate to high upstream Lorenz factor, and might alter the details of 
the resultant spectrum. 



6 Conclusions 

In this paper we consider sub-photospheric emission by relativistic radiation 
mediated shocks in the prompt phase of GRBs. The main conclusions are: 

1. Dissipation via internal shocks is likely to occur over a range of radii 
that encompasses the photosphere. Slower shells should collide below the 
photosphere, at a moderate optical depth (unless the time separation be- 
tween consecutive ejections is much larger than the dynamical time of the 
engine), whereas shells ejected with sufficiently large Lorentz factor will col- 
lide above the photosphere. The fraction of the energy that dissipates below 
the photosphere depends on the burst parameters, and may significantly vary 
from burst to burst. Specific examples, for which measurements of the lu- 
minosity and estimates of the Lorentz factors are available were analyzed. 
In case of GRB 090902B we conclude that a major fraction of the explosion 
energy dissipated below the photosphere, in a region of optical depths < 300, 
via relativistic radiation mediated shocks. For GRB 080916C we find that 
collision of shells occurred predominantly above the photosphere, giving rise 
to formation of collisionless shocks. 

2. Shocks that form below the photosphere are mediated by Compton 
scattering of radiation produced inside the shock or advected by the upstream 
flow. Since the scale of the shock, a few Thomson lengths, is vastly larger 
than any kinetic scale involved, particle acceleration by the Fermi process 
is highly unlikely in such shocks. On the other hand, bulk Comptonization 
gives rise to a hard spectrum extending above a (thermal) Wien peak, up to 
the KN limit •y u m e c 2 (and perhaps beyond, e.g., Budnik et al. 2010), where 
7„ is the Lorentz factor of the upstream flow, measured in the shock frame. 
The mean photon energy in the immediate downstream is determined by the 
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number of photons advected into and/or produced inside the shock. For a 
cold upstream, a condition expected in shock breakout episodes, the Wien 
temperature is ~ m e c 2 (Budnik et al. 2010). However, under the conditions 
anticipated in the prompt phase of GRBs, the shock upstream is photon rich, 
in the sense that photon advection largely dominates over photon production, 
and the thermal peak is located at energies well below m e c 2 . We estimate 
temperatures in the range between a few and a few tens keV for typical GRB 
parameters. 

The existence of the hard spectral component is demonstrated in this work 
through test particle Monte Carlo simulations on given shock profiles, that 
were computed in an earlier work (Levinson and Bromberg 2008), neglecting 
pair production by MeV photons and the thermal spread of the electrons. 
The resultant spectra extend to the KN limit, which imply that pair pro- 
duction inside the shock is likely to be important in photon rich shocks. A 
full treatment requires a proper account of this process in the computation 
of the shock structure and the resultant spectrum. 
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Figure 1 The optical depth at the shock location calculated for different 
Lorentz factors and time separation of the ejected shells, in 3 Fermi GRBs. 
The curve lines depict constant St, taken between a minimal dynamic time 
scale to a maximal value taken to be the observed variability by Fermi BGO 
detector. The colored areas show the range of T and St relevant for each 
burst. 




-^f) (orange), 

presented on a 7 U /3 U vs. N (left panels) and P (right panels). Areas marked 
in gray correspond those parts in the transition region where the advection 
of photons from the upstream dominates over emission processes, and the 
temperatures are NR. The green line depicts constant values of Zopio(T), 
where values < imply that a high energy spectral tail may evolve in the 
shock transition. Below the black line shocks can not form since the Mach 
number < 1. 




Figure 3 Contour lines of equal ( c (light blue) and log 10 ( —^f- ) (orange), 

presented on a r y u f3 u vs. r]/r] c plane. The values here are calculated on the 
photosphere. Areas marked in gray correspond those parts in the transition 
region where the photons advection dominates over emission processes, and 
the temperatures are NR. The green line depicts constant values of Zopio(T), 
where values < imply that a high energy spectral tail may evolve in the 
shock transition. Below the black line shocks can not form since the Mach 
number < 1. 
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Figure 4 The angle averaged spectral energy distribution in the immediate 
downstream of a RRMS of width r = 3 and Lorentz factor 7„ = 2 (red) and 
7„ = 7 (blue), as measured in the shock frame. The downstream velocity in 
both cases is taken to be (3d = 1/3. The calculated energy is normalized to 
the initial injection energy, 7de£ n , with e' in = 100 ev in the fluid frame. 




Figure 5 Same as Fig EJ but for an upstream Lorentz factor 7„ = 2 and 
different values of r. 



